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may have some bearing on the formation of primordial black holes in Randall-Sundrum models. 
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1. Introduction 

The Randall-Sundrum (RS) braneworld model ^ postulates that our observable universe is 
a thin 4-dimensional hypersurface residing in 5-dimensional anti-de Sitter (AdS) space. Ordinary 
matter degrees of freedom are assumed to be confined to the brane, while gravitational degrees of 
freedom are allowed to propagate in the full 5-dimensional bulk. The warping of AdS space allows 
us to recover ordinary general relativity (GR) at distances greater than the curvature radius of the 
bulk £. Current laboratory tests of Newton's law constrain £ to be less than around 0.1 mm 

It is well known that the Friedmann equation governing the expansion of the brane universe 
differs from general relativity by a correction of order p/o, where p is the energy density of the 
brane matter and a > (TeV)'* is the brane tension. The magnitude of this correction defines the 
"high-energy" regime of braneworld cosmology as the era when p > a or equivalently Hi > 1, 
where H is the Hubble parameter. 

The equations of motion governing fluctuations of the model are found to differ from GR 
in two principal ways at early times: First, they acquire 0'{p /a) high-energy corrections similar 
to those found in the Friedmann equation. By themselves, such corrections are not difficult to 
deal with: they just modify the second-order ordinary differential equations (ODEs) governing 
perturbations in GR. But the second type of modification is more problematic: perturbations on the 
brane are also coupled to fluctuations of the 5-dimensional bulk geometry, which are collectively 
known as the "Kaluza Klein" (KK) degrees of freedom of the model. The KK modes are governed 
by master partial differential equations (PDEs) defined throughout the AdS bulk |^ ||]. The only 
known way of solving this system of equations on all scales simultaneously is by direct numerical 
solution. 

The purpose of the paper is to numerically solve for the behaviour of scalar perturbations in 
the radiation-dominated regime of braneworld cosmology. We use two different numerical codes 
recently developed in Refs. [^, ^, which gives us the ability to confirm the consistency of our 
numeric results via two independent algorithms. We are ultimately interested in finding the matter 
transfer function in the radiation era, and also determining the relative influence of KK and high- 
energy effects on the density perturbations. Heuristically, we may expect the KK modes to amplify 
high-energy/small-scale density perturbations. The reason is that we know that the gravitational 
force of attraction in the RS model is stronger than in GR on scales less than £ |^ This implies 
that modes with a physical wavelength smaller than I during horizon crossing will be amplified 
due to the KK enhancement of the gravitational force. However, this physical reasoning needs to 
be tested with numeric simulations. 

2. Scalar perturbations 

It has been shown in Refs. [^, ^] that scalar-type perturbations of the bulk geometry are gov- 
erned by a single gauge invariant master variable Q.. This bulk master variable satisfies the follow- 
ing wave equation 




(2.1) 
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and a boundary condition on the brane 

The density contrast on the brane A satisfies the wave equation 

^2 i^2^2 



0. 



(2.2) 
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A = 6c,2 - 1 - 8w + 3w^ 



B = 3ci-9w- 4. 



(2.3a) 
(2.3b) 



Here, we have defined the sound speed = 8p/8p and the equation of state w = p/p (where 
p and p are the energy density and pressure of the brane fluid, respectively), F is the entropy 
perturbation of the brane matter and Tj is the conformal time along the brane. In this paper we 



assume that the matter anisotropic stress vanishes. The above ODE, the bulk wave equation (gjj) 



and the boundary condition ( |2.2D comprise a closed set of equations for A and D.. Note that in the 
low energy universe, we can neglect ^{p^/ cj2) terms. If we also set = we obtain the standard 
4-dimensional dynamical equation for A; hence, we recover GR at low energies. 



3. Numeric analysis and discussion 

For the rest of the paper, we will restrict ourselves to the case of a radiation-dominated brane 
with w = 1/3. We define the "*" epoch as the moment in time when a mode with wavenumber 
k enters the Hubble horizon, k = H^a^. Another important era is the critical epoch, the transition 
between high and low energy regimes, when Hcl = 1 and the radiation density has its critical value 
Pc/o = Vl — 1. Generally speaking, we call modes with k > k^ "supercritical" and modes with 
k < kc "subcritical". The scale defined by the critical mode in today's universe corresponds, for 
i = 0.1mm, to a scale of ~ 10 astronomical units (AU), which is incredibly tiny by cosmological 
standards. 

In Fig. H, we plot the predictions of GR, the 4-dimensional effective theory (where all &{p/o) 
corrections to GR are retained, but the bulk effects are removed by artificially setting D. = 0), 
and the full 5-dimensional simulations for the behaviour of the curvature perturbation on uniform 
density slices and the density contrast A for a supercritical mode. Since in any given model we 
expect the primordial value of the curvature perturbation to be fixed by inflation, it makes physical 
sense to normalize the waveforms from each theory such that i^^f, « i^qr ~ 1 for a ^ a*. 

When this is enforced we see that the effective theory predicts a larger final amplitude for the 
density perturbation than GR. Furthermore, the final amplitude in the 5-dimensional simulation is 
larger than both of the other theories. From this we can infer that, as we expected, both ^(p/a) 
and KK effects induce enhancement in the amplitude of perturbations. 

As in Fig. [T], let the final amplitudes of the density perturbation with wavenumber k be %o{k), 
'(Cff(^) and 'i^GRik) for the 5-dimensional, effective and GR theories, respectively. Then, we define 
enhancement factors as 
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Figure 1: A comparison of the behaviour of the curvature perturbation ^ (left) and the density perturbation 
A (right) in the full 5-dimensional theory including KK contributions (5D), the effective 4-dimension theory 
including G{p j a) corrections (eff), and ordinary general relativity (GR). The waveforms for each theory 
are normalized such that ^ = 1 on superhorizon scales. 



It follows that ^cff(^) represents the iff{p /o) enhancement to the density perturbation, £2s{k) gives 
the magnification due to KK modes, while ^s^ik) gives the total 5-dimensional amplification over 
the GR case. These enhancement factors are shown in the left panel of Fig. ^ We can see that 
they all increase as the scale is decreased, and that they all approach unity for k^O, which means 
we recover general relativity on large scales. For all wavenumbers we see £2^ff > £2g > 1 , which 
implies that the amplitude magnification due to the ^{p/a) corrections is always larger than that 
due to the KK modes. Interestingly, the .^-factors appear to approach asymptotically constant 
values for large k. 

Now we consider a transfer function T{k) that will tell us how the initial spectrum of curvature 
perturbations maps onto the spectrum of density perturbations at some low energy epoch 
within the radiation era characterized by the conformal time T] > T]c- It is customary to normalize 
transfer functions such that T(k;r\) 1, which leads us to the following definition 





T{k;n] 



iiif 



Here, Q"' is the primordial value of the curvature perturbation and ^k{'^ ) is the maximum amplitude 
of the density perturbation in the epoch of interest. We know that we recover the GR result in 
the extreme small scale limit, which gives the transfer function the correct normalization. In the 
righthand panel of Fig. ^ we show the transfer functions derived from GR, the effective theory and 
the 5-dimensional simulations. As expected, the T{k;r]) for each formulation match one another 
on subcritical scales k < k^. However, on supercritical scales we have T^o > T^n > 7gr. Our results 
are robust against modifications of the initial data for simulations. 

The amplitude enhancement of perturbations is important on comoving scales < 10 AU, which 
are far too small to be relevant to present-day/cosmic microwave background measurements of the 
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Figure 2: Density perturbation enhancement factors (left) and transfer functions (right) from simulations, 
effective theory, and general relativity. The transfer functions in the right panel are evaluated at a given 
subcritical epoch in the radiation dominated era. 



matter power spectrum. However, it may have an important bearing on the formation of compact 
objects such as primordial black holes and boson stars at very high energies, i.e. the greater grav- 
itational force of attraction in the early universe will create more of these objects than in GR. We 
know that the abundance of primordial black holes can be constrained by big bang nucleosythe- 
sis and observations of high-energy cosmic rays, so it would be interesting to see if the kind of 
enhancement of density perturbations predicted in this paper can be used to derive new limits on 
Randall-Sundrum cosmology. 
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